Emergence of long distance pair coherence through incoherent 
local environmental coupling 
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We demonstrate that quantum coherence can be generated by the interplay of coupling to an 
incoherent environment and kinetic processes. This joint effect even occurs in a repulsively in- 
teracting fermionic system initially prepared in an incoherent Mott insulating state. In this case, 
coupling a dissipative noise field to the local spin density produces coherent pairs of fermions. The 
generated pair coherence extends over long distances as typically seen in Bose-Einstein condensates. 
This conceptually surprising approach provides a novel path towards a better control of quantum 
many-body correlations. 
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In recent years, various experimental methods have 
been developed to dynamically generate non-trivial cor- 
relations in quantum materials. On the one hand, exter- 
nal electromagnetic fields have been used to photo-induce 
phase transitions in solid state materials pQ. For exam- 
ple, spin density wave order was induced in the normal 
state of a pnictide compound using femtosecond opti- 
cal pulses P]. A Josephson plasmon, typically present 
in a superconducting state, has even been triggered in 
a non-superconducting striped-order cuprate by the ap- 
plication of mid-infrared femtosecond pulses [3] . On the 
other hand, environmental tailoring [4J has been used to 
prepare highly entangled states such as a Bell state of 
two ions [5] or a Tonks-Giradeau-like state in a molecu- 
lar quantum gas [6]. In these examples, the realization 
of complex states relies on the same principle as optical 
pumping whereby atoms are prepared in so-called dark 
states immune to environmental coupling. 

We report on a complementary mechanism where dy- 
namical generation of coherence is achieved through the 
combined effect of a simple local dissipative coupling and 
kinetic processes. To examplify the inner workings of this 
mechanism, we consider an ultracold fermionic gas in an 
optical lattice subjected to local spin-polarization mea- 
surements carried out by phase-contrast imaging [7] or 
spatial and temporal light field fluctuations. The dissi- 
pative coupling heats up the system and destroys single- 
particle correlations. At the same time the number of lo- 
cal pairs, which are immune to the dissipative coupling, 
increase due to kinetic hopping. Unexpectedly, these lo- 
cal pairs then act as a source for the generation of pair 
correlations over longer distances. The produced cor- 
relations are long lived and reminiscent to those of the 
celebrated ^-pairing state [HI [H] > a condensate of bound 
on-site pairs of momentum k — — (with a the lattice 
spacing). Moreover, the appearance of a sharp feature 
in the pair momentum distribution, as shown in Fig. [T] 
serves as a signature for the formation of long distance 
coherence. In cold atom experiments, such pair momen- 



tum distributions can be observed by the projection of 
the local pairs onto molecules (TUJ [TT] . As our proposal 
relies both on dissipation and kinetic processes, it is con- 
ceptually very different from previous approaches where 
the ?7-pairing state was stabilized through either adia- 
batic state preparation [91 [T2] , or the imprint of phase 
coherence between neighboring sites by a tailored envi- 
ronment [11 113j. 

The system under consideration here is made of re- 
pulsively interacting fermions on a <i-dimensional lattice 
of volume V and lattice constant a. We describe this 
many-body system by the Hubbard model 

H = - J Yl (4,<A V + h-c) + «r,tnr4 . 

(r,r'),er r 

where c\ a is the creation operator for a fermion with 
spin a =t , I and site index r, n r , CT = cj CT c r a is the 
density operator, J > is the hopping coefficient, U the 
interaction strength, and (r, r') indicates that the sum 
is done over nearest-neighbors. This Hamiltonian is one 
of the simplest models capturing the interplay between 
the kinetic and interaction energies, and can be used, 
for example, to understand the metal to Mott insulator 
transition. A particularly clean realization of this model 
is achieved using ultracold fermionic gases confined to 
optical lattices [H] . 

In the present work, we assume the system to be ini- 
tially prepared in a stationary state of this Hamiltonian, 
typically a Mott insulator as realized in Refs. [13 fTo] . 
We study the system dynamics after coupling a dissipa- 
tive environment to the local spin densities such that 

±p(t) = -^[H,p(t)]+V[p{t)}, (1) 

with 



V [p{t)\ = ( n s ^ r p h 8>r - 
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FIG. 1: Time evolution of the momentum distribution of lo- 
cal pairs. A chain of 36 lattice sites is prepared at t = in a 
perfect Mott insulating state where pair correlations are ab- 
sent. A fast build-up in occupation of the momenta except 
close to k — — takes place at short times (plotted on a lin- 
ear time scale). Then, over time, all momenta except k = — 
become homogeneously occupied, signaling the generation of 
the coherence of pairs over longer distances (plotted on a log- 
arithmic scale). The evolution is obtained using the effective 

diffusion equation Q with ^ = 1.5 and 70 = j^f- 



The quantum jump operators n s r = r — h\, tT measure 
the local spin polarization. The dissipative coupling T> 
can be realized, for example, by a light field whose fre- 
quency is chosen in between the transitions of the two 
fermionic states as used for phase contrast imaging [TJ. 
This light field can either be used to probe the spin den- 
sity locally, e.g. in combination with an independent ad- 
dressing of each site p~7j [18] or to create "magnetic field" 
noise, e.g. through the realization of spatially-disordered 
and time-decorrelated white noise patterns. 

This dissipative mechanism leads to an expo- 
nential decay of single particle correlations as 
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In contrast, pair correlations 

c a = v J2 r , remain unchanged 

under the action of the dissipator, T>, as doublons 
(doubly occupied sites) and holes (empty sites), which 
have no net polarization, belong to the dissipation-free 
subspace. In particular, the ^-paired state, gener- 
ated through the repeated application of the operator 

ff* = J2 r e%7V r ^rt^4 011 ^ e vacuum J is P ar1, °f this sub- 
space; here 7r = (71", . . . , 7r). As the 77-paired state is an 
eigenstate of the Hamiltonian, it is immune against the 
action of both the unitary and dissipative operators (see 
Eq. 0). Interestingly, even if the initial state does not 
overlap with the 77-paired state, correlations may emerge 
due to kinetic processes. 

Irrespective of the coupling strength and the properties 
of the Hamiltonian, at sufficiently large times, Tt ^> 1, 
the dissipation free subspace is reached. This subspace 
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FIG. 2: Left: Example of the effective creation and diffusion 
of pair correlations. Right: Within the adiabatic elimination 
method, the evolution is based on the effective coupling of 
two states (lower and upper state) of the decoherence free 
subspace via a virtual excitation (center). The virtual state 
is reached by the hopping process and can decay with V and 
dephase due to interaction U. Here this process is examplified 
for a state with no pair correlations (lower state) connected 
to a state containing pair correlations (upper state) through 
the creation process (box) presented on the left panel. 



is highly degenerate with respect to the dissipator T> 
and the Hamiltonian can lift this degeneracy. To under- 
stand the dynamics, we perform adiabatic elimination 
(see Appendix) revealing how hopping-induced virtual 
excitations, around the dissipation-free subspace, affect 
the evolution of the system (cf. Fig. |2j. The effective 
coupling via the virtual excitations depends on whether 
the interaction energy is changed during the process and 
takes the form 



7o 



8J^ 
h 2 T 



and 



lu 



8J 2 r 
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Using this perturbative approach, the equations de- 
scribing the evolution of staggered pair correlations, 
i.e. Cd = — e l7r d Cd H — ip-, for times larger than ^ 
are cast into a system of coupled diffusion equations 
(cf. Fig. D left panel): 



,|d-d'|=i 



d' jd (t) (c d ,(t)-c d (t) 



(3) 



The diffusion constant depends on the coupling to the 
different virtual excitations weighted by their probability 
to occur, 



A d ,, d (f) = D(t) =70 7 + °o(t) + lu 



C (t) 



for |d| and |d'| 5= 0, while A A , t0 {t) = A QA (t) = \. 
For the sake of concreteness, we assumed above that the 
system was half filled and translationally invariant. Gen- 
eralizations are straightforward and do not lead to quali- 
tative changes. Moreover, we provide numerical evidence 
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that this diffusive behavior is not restricted to the domain 
r >• i , but is valid even at weak coupling T < jr. How- 
ever, in the latter case, the diffusion constant deviates 
from the perturbative results. 
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FIG. 3: Cd(t) as a function of time t as described by the 
diffusion equation Q . The chain is prepared in a perfect Mott 
insulator with L = 36 sites and evolves with ^ = 1.5. Arrows 
mark the t — > oo limit. The double occupancy, Cd=o, and 
nearest-neighbor pair correlation, Cd=i, raise quickly and feed 
the delayed increase of correlations at large distances. Inset: 
evolution of the diffusion constant as a function of time. D(t) 
becomes time-independent as the double occupancy, Cd=o, 
saturates. 

We illustrate the creation of correlations using, as an 
example, a system initially in a Mott insulating state, 
i.e. Cd(0) = for all d. In Fig. [3] we depict the dynam- 
ics triggered by the action of the Hamiltonian, H, and 
the dissipator, T>, on a one-dimensional system. First, 
double occupancy and short range pair correlations rise 
on the time scale — . Then, following this initial build- 
up, the double occupancy and the nearest-neighbor pair 
correlation act as sources for the propagation of pair 
correlations over longer distances. Within the pertur- 
batively derived Eq. ([3]), one expects the propagation 
of the staggered pair correlations to be described by a 
normal diffusion process. The evolution of the system 
is constrained by the constant of motion F = C^ =n / a , 
where Ck = y J^d e~ Cd is the momentum distri- 
bution of the local pairs (Fig.[T]). The constant of motion 
is proportional to the number of 77-pairs, (fj'f]}, initially 
present in the system. For an atomic Mott insulator, 
(v v) = 0- Thus, the pair correlations asymptotic val- 
ues are C d {t -> 00) = - e l7r d ^ for |d| ^ while 
C (t -> 00) = \ (1 - i). The "overall" sign of C* d ^o, 
and thus the phase slip between d = and d = 1 de- 
pends on the initial state and would differ if F > \. We 
note that the saturation of the double occupancy implies 
that the diffusion constant becomes time- independent: 
D(t ->■ 00) ~ :m± ^ L (see inset in Fig. [3}. To study this 
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FIG. 4: Cd(t), the staggered pair correlations, are shown 
as a function of distance, d, for three different times in a 
chain of L = 36. We compare the solution of the diffusion 
equation Q with time-dependent DMRG simulations. For 
both cases the initial conditions are taken from ground state 
DMRG calculations at U — 12 J. We use |p = 1.5 and in the 
time-dependent DMRG U = 12J. Inset: Symbols represent 
the variance of the pair correlation distribution versus time. 
Lines are linear fits for 10 < 7ot < 20. 

diffusive spreading, we plot in Fig. [4] the staggered corre- 
lations at different times versus distance. In terms of Cd, 
the initial Mott insulator is characterized by a peak at 
d = 0. As shown on Fig. |4j the interplay of dissipation 
and hopping gradually transforms this peak into a broad 
gaussian distribution. 

To check the validity of our perturbative results against 
unbiased methods, we use time-dependent DMRG to 
solve the master equation ([I]) stochastically [20] . One can 
appreciate in Fig. [4] that the pair correlations obtained 
from DMRG are in good agreement with the pertubative 
results. Deviations away from the expected gaussian dis- 
tribution mostly occur at short distances. This discrep- 
ancy is partially attributed to the decoupling of density 
correlations applied in the derivation of Eq. ([3| . 

The diffusive propagation is best characterized by the 
variance 

d^ = ]Td 2 Cd(t)/£<7d(o). 

d d 

Generally, for a hypercubic lattice with connectivity z, 
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terns. In contrast to correlations realized in cooled con- 
densed matter systems, typically sensitive to tempera- 
ture, the non-equilibrium mechanism presented above is 
immune against thermal fluctuations. This conceptually 
surprising approach provides a new route towards a bet- 
ter control of quantum many-body correlations. 
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FIG. 5: The variance of the distribution of the staggered 
pair correlations Cd(t) with distance for different dissipative 
strength ^- at constant U = 12J in a chain of length 36. 
The linear rise of the variance at larger times is in accordance 
with a diffusive formation of the pair correlations. Symbols 
are taken from DMRG simulations, straight lines are linear 
fits d 2 = 2_D e ff t + const, for values with d 2 > 2.5. Inset: D c ff 
as a function of dissipative coupling at U = 12 J (symbols). 
The straight line corresponds to D(t — > oo) derived within 
perturbation theory and agrees well with the numerical re- 
sults for large coupling strengths. 



the variance of a diffusive process obeys 



d_ 
dt 



d 2 (t) = zD(t) 



Within the perturbative treatment, at t 3> where D(t) 
becomes constant, the variance should rise linearly with 
time. The variance, shown in the inset of Fig. [4] confirms 
this statement. Small deviations from the linear behavior 
are consistent with the time-dependence of D{t). 

Remarkably, we find the diffusive description of the 
propagation of pair correlations to remain valid down to 
the weakly dissipative regime. However, in this regime, 
the diffusion constants need to be phcnomcnologically 
determined. In Fig. [5j we study a strongly interacting 
system U = 12J with couplings from T = 4jr down to 
r = 0.25^. In all cases, a normal diffusive regime is 
entered after t ~ jj. As shown in the inset of Fig. [HJ 
the effective diffusion constants, in the strong dissipative 
coupling limit, agree nicely with our perturbative pre- 
dictions. As expected, with decreasing T, the deviation 
from the analytic predictions increases. Nevertheless, the 
correct qualitative behavior is predicted for T > ^: the 
effective diffusion constant increases with i. For T < j- 
our simulations suggest a saturation of the diffusion con- 
stant to roughly 0.7 j-. 

In summary, we demonstrated here that the combined 
action of incoherent local environmental coupling and ki- 
netic processes can result in the emergence of long dis- 
tance pair correlations in repulsive fermionic lattice sys- 



Appendix: Adiabatic elimination 

Locally, the dissipation free subspace with respect 
to T> is spanned by the diagonal operators {p r q — 

\0)(0\,Pr,l = imi,Pr,t = lt)<tj,Pr,n = \U)(U\h 

and the off-diagonal operators {d r = |0)(tI|,<4 = 
|t!)(0|} annihilating or creating a pair at site r 
such that dl — cj^cj^. In addition, the first ex- 
cited subpace can be defined via the basis elements 
{Pr,oc r ,<T , 4,<rPr,o , Pr,4.f4,a > c r ,aPr,n}- These operators 
form a diagonal basis for — [p r ,ti; ■] +2? [.] with eigen- 



values X r 



1. 



c l 2 ' 2 ' 2 ^ 4 ft ' 2 ' H > ' 

We associate a projector V" to each subspace associated 
with a certain eigenvalue of the excited subspace, while 
projects onto the dissipation free subspace. Via adi- 
abatic elimination of the excited subspace (see e.g. [19]). 
one derives effective equations of motion for the basis el- 
ements within the dissipation free subspace, for example 



E 



K r y,V?V?, 



K i d 



dt~* ^ h? 

r' , |r— r' | — 1 
a.,a' — 1... . ,4 



with A r 



and K r 



r' . |r— r'l — 1 



'V 1 r 



(A.4) 



Pr,0 



c 



h.c. Eq. (A.4 1 is used to derive the 



equation of motion for d\d r i with |r — r'| > 1. A similar 
procedure is used to find the equations for |r — r'| < 1. In 
order to construct the closed set of equations of motion 
for the pair correlators, Eq. ([3|, we decouple (p r ,nOr') as 
(Pr,n)(0 T '} where n = {0,t,!,t!} an d O r > is an arbitrary 
local operator on site r' ^ r. 
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